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In the Israel-Stewart's theory of 2nd order dissipative hydrodynamics, we have simulated <f> pro- 
duction in Au+Au and Cu+Cu collisions at -^^^=200 GeV. Evolution of QGP fluid with viscos- 
ity over the entropy ratio ?7/s=0.25, thermalised at Ti=0.2 fm, with initial energy density Ei=5.l 
GeV/fm 3 explains the experimental data on 4> multiplicity, integrated V2, mean pr, pr spectra and 
elliptic flow in central and mid-central Au+Au collisions. g/s=0.25 is also consistent with centrality 
dependence of <j> pr spectra in Cu+Cu collisions. The central energy density in Cu+Cu collisions is 
£i =3.48 GeV/fm 3 . 
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I. INTRODUCTION 

Experiments in Au+Au collisions at RHIC [l], 0, d, 
4], produced convincing evidences that in non-central 
Au+Au collisions, a hot, dense, strongly interacting, col- 
lective QCD matter is created. Whether the matter 
can be characterized as the lattice QCD [5, 6] predicted 
Quark-Gluon-Plasma (QGP) or not, is still a question of 
debate. For long, strangeness enhancement is considered 
as a signature of QGP formation Q. In QGP environ- 
ment, gg — > ss is abundant. If not annihilated before 
hadronisation, early produced strange and anti-strange 
quarks will coalesce in to strange hadrons and compared 
to elementary pp collisions, strange particle production 
will be enhanced. Recently, STAR collaboration pub- 
lished their measurements of 4>(ss) mesons in Au+Au 
[1, [9[ and in Cu+Cu [HI collisions. Both in Au+Au and 
Cu+Cu collisions, compared to pp collisions, (f> meson 
production is enhanced. However, it is uncertain whether 
or not the enhancement is due to increased production in 
QGP or due to canonical suppression of strangeness in pp 
collisions. STAR measurements of (j> mesons in Au+Au 
collisions found to be compatible with a model based on 
recombination of thermal s quarks [TTI | , strengthening the 
belief that in Au+Au collisions, a robust thermal parton 
source is created. 

Relativistic hydrodynamics provides a convenient tool 
to analyse Au+Au collision data. It is assumed that in 
the collision a fireball is produced. Constituents of the 
fireball collide frequently to establish local thermal equi- 
librium sufficiently fast and after a certain time r< 5 hydro- 
dynamics become applicable. If the macroscopic proper- 
ties of the fluid e.g. energy density, pressure, velocity 
etc. are known at the equilibration time Tj, the relativis- 
tic hydrodynamic equations can be solved to give the 
space-time evolution of the fireball till a given freeze-out 
condition such that interactions between the constituents 
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are too weak to continue the evolution. Using suitable 
algorithm (e.g. Cooper-Frye) information at the freeze- 
out can be converted into particle spectra and can be 
directly compared with experimental data. Thus, hydro- 
dynamics, in an indirect way, can characterize the initial 
condition of the medium produced in heavy ion collisions. 
Hydrodynamics equations are closed only with an equa- 
tion of state (EOS) and one can investigate the possibility 
of phase transition in the medium. A host of experimen- 
tal data produced in Au+Au collisions at RHIC, at cm. 
energy -^5=200 GeV, have been successfully analysed us- 
ing ideal hydrodynamics [12| , with an equation of state 
with 1st order confinement-deconfinement phase transi- 
tion. Multiplicity, mean px, py-spectra, elliptic flow etc. 
of identified particles, are well explained in the ideal hy- 
drodynamic model with QGP as the initial state. Ideal 
hydrodynamics analysis of the RHIC data indicate that 
in central Au+Au collisions, at the equilibration time 
T, w 0.6 fm, central energy density of the QGP fluid is 
£i «30 GeV I fm~ 3 12]. It may be mentioned that ideal 
hydrodynamics description of data are not unblemished. 
Pt spectra or the elliptic flow are explained only up to 
transverse momenta px ~ l.5GeV. At higher px descrip- 
tion deteriorates. Also ideal hydrodynamic description to 
data gets poorer in peripheral collisions. 

However, estimate of initial condition of the fluid can 
not be creditable unless dissipative effects are accounted 
for. Unlike in ideal fluid evolution, where initial and 
final state entropy remains the same, entropy is gener- 
ated in viscous evolution. Consequently, to produce a 
fixed final state entropy, viscous fluid require less ini- 
tial energy density than an ideal fluid. QGP viscosity is 
quite uncertain. Theoretical estimate cover a wide range, 
rj/s w 0-1. String theory based models (ADS/CFT) give 
a lower bound on viscosity of any matter rj/s > 1/An 
[l3|]. In a perturbative QCD, Arnold et al [l4[ esti- 
mated r]/s ~ 1. In a SU(3) gauge theory, Meyer [l5[ 
gave the upper bound T]/s <1.0, and his best estimate is 
r?/s=0. 134(33) at T = 1.165T C . At RHIC region, Naka- 
mura and Sakai [l6| estimated the viscosity of a hot gluon 
gas as r]/s=0. 1-0.4. Attempts have been made to es- 
timate QGP viscosity directly from experimental data. 
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Gavin and Abdel-Aziz [l7| proposed to measure viscos- 
ity from transverse momentum fluctuations. From the 
existing data on Au+Au collisions, they estimated that 
QGP viscosity as ?7/s=0.08-0.30. Experimental data on 
elliptic flow has also been used to estimate QGP viscos- 
ity. Elliptic flow scales with eccentricity. Departure form 
the scaling can be understood as due to off-equilibrium 
effect and utilised to estimate viscosity [l8[ as, rj/s—0. 11- 
0.19. Experimental observation that elliptic flow scales 
with transverse kinetic energy is also used to estimate 
QGP viscosity, tj/s ~ 0.09 ± 0.015 [19], a value close 
to the ADS/CFT bound. From heavy quark energy 
loss, PHENIX collaboration [2(| estimated QGP viscos- 
ity tj/s « 0.1-0.16. 

In recent years, considerable progress has been made in 
numerical implementation of dissipative hydr ody namics 
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the viscous hydrodynamic simulation of elliptic flow in 
Au+Au collisions, Luzum and Romatschke 34j obtained 
an upper bound to the ratio r\j s < 0.4-0.5. In (3f| Song 
and Heinz also argued that -q/s < 0.4 is a robust upper 
bound of QGP viscosity. In a recent paper (36|, we have 
estimated QGP viscosity as r\js w0.25. It was shown 
that 4> mean px is sensitive to QGP viscosity. STAR 
data [1, Q on centrality dependence of mean px in 
Au+Au collisions definitely reject ideal fluid or fluid with 
viscosity r\j s < 0.08-0.16. Data are explained only with 
??/s=0.25. Initial central energy density of the fluid is 
£i=5.1 GeV, much less than that the estimated value 
~ 30 GeV j fm in ideal hydrodynamics. Evolution of 
viscous fluid (?7/s=0.25) also explains the centrality de- 
pendence of <f> multiplicity, integrated i>2, Pt spectra up 
to w 3 GeV. 

Purpose of the present paper is to show that viscos- 
ity over entropy ratio ?7/s=0.25, is consistent with the 
STAR measurements of <p elliptic flow in central and mid- 
central Au+Au collisions. It is also consistent with the 
recent STAR data on </> production in Cu+Cu collision 
at 200 GeV. The central energy density however, is less 
in Cu+Cu collisions, £^=3.48 GeV/ fm 3 . The paper is 
organised as follows: in sectionfllj we briefly describe the 
hydrodynamical equations used to compute the evolution 
of ideal and viscous fluid. We have used a lattice moti- 
vated equation of state. Construction of the EOS is also 
discussed in section [TTJ Simulation results are discussed 
in section IIII1 Summary and conclusions are given in 
section ITVl 



II. HYDRODYNAMICAL EQUATIONS, 
EQUATION OF STATE AND INITIAL 
CONDITIONS 

A. Hydrodynamical equations 

In the Israel-Stewart's theory of 2nd order dissipative 
hydrodynamics, space-time evolution of the fluid is ob- 
tained by solving, 
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FIG. 1: Viscous fluid (r?/s=0.08) simulation for temporal 
evolution of momentum anisotropy in b=7 fm Au+Au col- 
lision at RHIC. The solid line is the simulation result from 
VISH2+1 [3;| and the dashed line is the simulation result 
from AZHYDRO-KOLKATA. Initial condition of the fluid is 
very similar in both the simulations. 
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EqQ] is the conservation equation for the energy- 
momentum tensor, T^ v = (e + p^u^u 11 — pg^ v + tt^ , 
e, p and u being the energy density, pressure and fluid 
velocity respectively. tt^ v is the shear stress tensor (we 
have neglected bulk viscosity and heat conduction). EqfJ] 
is the relaxation equation for the shear stress tensor tt^ . 
In EqfJl D = u^d^ is the convective time derivative, 
v <m^> = !(v*V + VV) - \{d.u){g^ - u^u") is a 
symmetric traceless tensor, rj is the shear viscosity and 
tv is the relaxation time. It may be mentioned that in 
a conformally symmetric fluid relaxation equation can 
contain additional terms [33| . 

Assuming boost-invariance, EqsQ] and [2] are solved in 
(t = \/t 2 — z 2 ,x,y,T) a = 2 m fr?) coordinates, with a 
code " 'AZHYDRO-KOLKATA'", 2 developed at the Cy- 
clotron Centre, Kolkata. Details of the code can be 
found in [221. To show that AZHYDRO-KOLKATA com- 
putes the evolution correctly, in Fig[TJ we have com- 
pared the temporal evolution of momentum anisotropy 

£p ^rp xx j^rpyy y 

Song and Heinz [3_ 
both the simulations. Within 10% or less, AZHYDRO- 
KOLKATA simulation reproduces Song and Heinz 's [33[ 
result for temporal evolution of momentum anisotropy 
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FIG. 2: Black circles are lattice simulation [y] for entropy 
density. The black line is the parametric representation to 
the lattice simulations. In the inset, the solid and dashed 
lines are the squared speed of sound in lattice based EOS and 
in an EOS incorporating 1st order transition [121 ] . 
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FIG. 3: (color online) The black and red solid lines are initial 
central energy density as a function of number of participants 
in Au+Au and Cu+Cu collisions. The dashed line shows 
the critical energy density for the confinement-deconfinement 
cross-over transition. 



B. Equation of state 

One of the most important inputs of a hydrodynamic 
model is the equation of state (EOS). Through this input 
macroscopic hydrodynamic models make contact with 
the microscopic world. Most of the hydrodynamical cal- 
culations are performed with EOS with a 1st order phase 
transition. Huovinen [371 ] reported an 'ideal' hydrody- 
namic simulation with 2nd order phase transition. He 
concluded that the experimental data (e.g. elliptic flow 
of proton or antiproton) are better explained with EOS 
with 1st order phase transition than with EOS with 2nd 
order phase transition. However, lattice simulations [6] 
indicate that confinement to deconfinement transition is 
a cross over, rather than a 1st or 2nd order phase transi- 
tion. It is then essential that hydrodynamic simulations 
are done with EOS with cross-over transition rather than 
with EOS with 1st or 2nd order transition. In FigfJl a 
recent lattice simulation [6( for the entropy density is 
shown. The solid line in Figf2] is a parameterisation of 
the entropy density. 

-^ = a+[p + 1 T}[l + tanh^^}, (3) 

From the parametric form of the entropy density, pres- 
sure and energy density can be obtained using the ther- 
modynamic relations, 

p(T) = [ s(T)ds (4) 
Jo 

e{T) = Ts-p. (5) 

Generally, in hydrodynamic simulations, hadronic 
phase is approximated by a (non-interacting) resonance 



hadron gas comprising all the resonances below 2-3 GeV. 
As the lattice simulation cover a wide temperature range 
below the cross over temperature, T co = 196(3) MeV, we 
choose to use the lattice based EOS (EqOO both in the 
QGP and in the hadronic phase. The idea is to expose 
the lattice simulation of EOS to experimental scrutiny. 
Indeed, lattice simulations are vague about the nature of 
the confined (T < T co ) phase. The confined phase is cer- 
tainly unlike hadronic resonance gas. The trace anomaly 
(e — 3p) /T 4 in the temperature range 140-200 MeV is ap- 
proximately 30% less than that of a hadronic resonance 
gas [6j. One also note that at low temperature, effective 
degrees of freedom in the confined phase is g^ sa2. In 
contrast, in hadronic resonance gas, gu > 40. In the in- 
set of Figf^l the squared speed of sound (cf ~ p/e) in the 
lattice based EOS is compared with c 2 s in an EOS with 
1st order phase transition which model the quark 
phase with bag model, and the hadronic phase by the 
hadronic resonance gas. In 1st order EOS, fall sharply 
near the critical temperature. The fall is smoothened out 
in cross over transition. Lattice based EOS is also softer. 



C. Initial conditions 

Solution of partial differential equations (Eqs lll2p re- 
quires initial conditions, e.g. transverse profile of the 
energy density (e(x,y)), fluid velocity (v x (x,y),v y (x,y)) 
and shear stress tensor (^"(x, y)) at the initial time Tj. 
One also need to specify the viscosity (rj) and the relax- 
ation time (ttt). A freeze-out prescription is also needed 
to convert the information about fluid energy density 
and velocity to particle spectra and compare with ex- 
periment. 
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FIG. 4: The top two panels show the energy density con- 
tours in 0-10% and 50-60% centrality Au+Au collisions. The 
contours are drawn at £=1.5 and 0.05 GeVj fm 3 . They corre- 
spond to confinement-deconfinement cross-over and freeze-out 
respectively. The bottom two panels show the energy density 
contours in Cu+Cu collisions. In 50-60% centrality Cu+Cu 
collisions, initially QGP is not produced. 



In [36j, we assumed that the fluid is thermalised at 
Ti=0.2 fm and the initial fluid velocity is zero, v x (x, y) = 
v y (x,y) — 0. Initial energy density was assumed to be 
distributed as [13| 



e(b,x,y) = e Q [0.75N part (h,x,y) + 0.257V co/; (b, x, y)}, 



(6) 

where b is the impact parameter of the collision. N part 
and N co u are the average participant and collision num- 
ber respectively. The shear stress tensor was initialised 
with boost-invariant value. For the relaxation time, we 
used the Boltzmann estimate 7V = Srj/Ap. The freeze-out 
was fixed at Tp=150 MeV. In Eqj6l £o is a parameter 
which does not depend on the impact parameter of the 
collision. Assuming that r\j s remain a constant through- 
out the evolution, for a set of values r)/s=0 (ideal fluid), 
0.08, 0.16 and 0.25, we fit e to reproduce STAR mea- 
surements of 4> multiplicity in 0-5% centrality Au+Au 
collisions. Centrality dependence of <fi multiplicity and 
<fi mean px are simultaneously explained only with vis- 
cosity over entropy ratio ?7/s=0.25. The corresponding 
central energy density, in b=0 Au+Au collision is e,;=5.1 
GeV/fm 3 . 

For Cu+Cu collisions, we only change the central en- 
ergy density, other parameters remain unchanged. The 
initial energy density in Cu+Cu collisions is obtained 
by fitting <j> multiplicity in central 0-10% Cu+Cu col- 
lisions. The fitted value corresponds to central energy 
density £^=3.48 GeV/fm 3 in b=0 Cu+Cu collisions. 
As it will be shown below, QGP fluid, initialised with 
£;=3.48 GeV/fm 3 reproduces most of the STAR mea- 



surements on cj) mesons in Cu+Cu collisions. It is inter- 
esting to note that e""/e™ » A x Ju /A\u . Apparently, 
initial energy density scales with nuclear radius. With 
boost-invariance, the initial system is a cylinder of radius 
R = tqA 1 / 3 , infinitely extended in the rapidity direction. 
Phenomenological relation ef u /ef u w A x Ju /A\u indicate 
that at the freeze-out also the system can be approxi- 
mated again by a cylinder with radius R' oc R. 

Before we compare hydrodynamic simulation for <f> pro- 
duction with experiment, it is interesting to compare ini- 
tial central energy density in different centrality rages of 
Au+Au and Cu+Cu collisions. In Fig[31 we have com- 
pared the initial central energy density in Au+Au and 
in Cu+Cu collisions as a function of participant num- 
ber. The black and red lines are for Au+Au and Cu+Cu 
collisions respectively. In the region where they overlap, 
initial central energy density in Au+Au and Cu+Cu col- 
lisions are similar. One then expects that <f> spectra in 
mid central Au+Au and central Cu+Cu collisions will 
be similar. The expectation is fulfilled in STAR experi- 
ment (e.g. 4> p T spectra in 40-50% Au+Au and 10-20% 
Cu+Cu collisions are nearly identical). In Fig[3J the 
dashed line is the energy density (e co wl.5 GeV/fm 3 ) 
for the confinement-deconfinement cross-over. In most 
of the collisions, in the central region, initially the fluid 
is produced in the deconfined state. However, energy 
density has a distribution, fluid in the central region is 
at higher density than the fluid at periphery. Thus in 
mid-central collisions, only a small portion of the fluid 
will be in the deconfined phase. In Fig in four pan- 
els, we have shown the contours of initial energy density 
in 0-10% and 50-60% centrality Au+Au and in Cu+Cu 
collisions. Contours are drawn at e co =1.5 GeV/fm 3 and 
£/ o =0.05 GeV/fm 3 , corresponding to cross-over energy 
density and freeze-out. In 0-10% centrality Au+Au and 
Cu+Cu collisions, initially, fluid in the central region is 
in QGP state. But fraction of fluid in QGP state is larger 
in Au+Au than in Cu+Cu collisions. One can immedi- 
ately say that hard probe signature of QGP formation 
will be less prominent in 0-10% centrality Cu+Cu col- 
lisions than in 0-10% centrality Au+Au collisions. For 
example, one can conjecture that in 0-10% Cu+Cu colli- 
sions J/t/>'s will be less suppressed than in a 0-10% cen- 
trality Au+Au collision. Experiments do vindicate the 
conjecture [H, HH, E3| • 



III. RESULTS 

A. centrality dependence of <j> multiplicity, mean pr 
and integrated V2 in Au+Au and Cu+Cu collisions 

In (36|, we have shown the viscous (?y/s=0.25) hydro- 
dynamics fit to the STAR data on <j) multiplicity, mean 
Pt and integrated flow in Au+Au collisions. For com- 
pleteness purpose, here also, we show the fits along with 
the fit obtained to Cu+Cu data. In Figs[5^ and[Sj3, the 
STAR measurements for the centrality dependence of (j> 
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FIG. 5: (a) Filled circles are STAR data on the centrality 
dependence of <j> meson multiplicity. The black line is viscous 
(tj/s=0.25) hydrodynamic fit to the data. The initial time 
71=0.2 fm, initial central energy density ei—5.1 GeV / fm 3 , 
freeze-out temperature is Tf = 150 MeV. (b) same as in (a) but 
for Cu+Cu collisions. The initial energy density is £;=3.48 
GeV/fm 3 . 



multiplicity {dN/dy) in Au+Au @, ^ and Cu+Cu [To| 
collisions are shown. In the region where N par t over- 
lap, <j) multiplicity is nearly identical in Au+Au and in 
Cu+Cu collisions. In Figj5^ and ,b the solid lines are 
the hydrodynamic predictions for </> multiplicity. Evo- 
lution of viscous (?7/s=0.25) QGP fluid thermalised at 
Tj=0.2 fm and initialised with central energy density 5.1 
GeV/ fm 3 in Au+Au collisions and 3.48 GeV/fm 3 in 
Cu+Cu collisions, reproduces the data in all the cen- 
trality ranges of collisions. As indicated above, we have 
used only the most central collision data (0-5% in case 
of Au+Au collisions and 0-10% in case of Cu+Cu col- 
lisions) to fix the initial energy density. Glauber model 
initial condition (Eq[5]) correctly incorporate the central- 
ity dependence and <p multiplicity is reproduced in all the 
centrality ranges of collisions. 

In FigJ5^, and b, we have shown the STAR measure- 
ments of<f> mean in Au+Au [1,0 and Cu+Cu [l(| col- 
lisions. Within the error bars, centrality dependence of <p 
mean px in Au+Au and in Cu+Cu collisions are nearly 
identical (though central value is consistently higher in 
Au+Au collisions). The black lines in Fig[5^,b are fit to 
the data in viscous hydrodynamics. In viscous hydrody- 
namics also, (f> mean px do not show any appreciable de- 
pendence on system size. STAR measurements of mean 
Pt in Au+Au collisions are nicely reproduced. <j) mean 
Pt in Cu+Cu are reproduced within 10% or less. 

Black lines in Fig[7£i and b, are the viscous hydrody- 
namics predictions for the centrality dependence of inte- 
grated V2 in Au+Au and Cu+Cu collisions. In Au+Au 
collisions, STAR measured integrated V2 in 0-5%, 10-40% 
and 40-80% centrality collisions [1, Filled circles in 
FigJT^ are the STAR measurements. Except for the very 
peripheral collision, hydrodynamic prediction is agree- 
ment with STAR data. In Cu+Cu collisions, integrated 
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FIG. 6: (a) Filled circles are STAR data Q on centrality de- 
pendence of mean pr of <fr mesons in Au+Au collisions. The 
black line is the fit obtained to the data in viscous hydrody- 
namics, (b) same as in (a) but for Cu+Cu collisions. 
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FIG. 7: (a) filled circles are STAR data on integrated elliptic 
flow in Au+Au collisions. The black line shows the centrality 
dependence of integrated V2- (b) centrality dependence of 
integrated vi in Cu+Cu collisions. 



i>2 is not measured yet. However, simulation results in- 
dicate that compared to Au+Au collisions, 4> meson in- 
tegrated flow is less in Cu+Cu collisions. The reason is 
understood. Initial eccentricity is small in Cu+Cu than 
in Au+Au collisions. Elliptic flow has size dependence, 
smaller the system, less is the flow. 



B. <j> pr-spectra in Au+Au and Cu+Cu collisions 

STAR measurements Q for <fi meson px spectra in 0- 
5%, 5-10%, 10-20%,20-30%, 30-40%, 40-50%, 50-60% and 
60-70% centrality Au+Au collisions are shown in Figj8^i. 
As noted by the STAR collaboration Q (f> py-spectra up 
to 30-40% centrality collisions are well fitted by an ex- 
ponential, indicating thermal production of <f> in central 
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FIG. 8: (color online) (a) STAR data Q tj> meson pr spectra 
in 0-5%, 0-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60% 
and 60-70% centrality Au+Au collisions. The black lines are 
4> spectra from evolution of viscous fluid, (b) STAR data 
on pt spectra of <f) in Cu+Cu collisions in 0-10%, 10-20%, 
20-30%, 30-40%, 40-50% and 50-60% Cu+Cu collisions. The 
black lines are viscous hydrodynamic predictions. 



collisions. Exponential fit worsen in more peripheral col- 
lisions. A Levy function (which has an exponential shape 
at low px and power law shape at large px) fits the pe- 
ripheral data. Apparently, in peripheral collisions non- 
thermal source contribute to <f> meson production. The 
black lines in FigHK are px spectra from evolution of vis- 
cous QGP fluid. Except for very peripheral collisions, 
data up to pr=3 GeV are well explained in viscous hy- 
drodynamics. At larger pt (not shown in FiglH]), viscous 
hydrodynamics under predict the py-spectra. At large 
Pt, other sources e.g. pQCD processes can contribute 
and hydrodynamic models may not be reliable. 

STAR collaboration recently published their measure- 
ments for (f> py-spectra in Cu+Cu collisions [1(3]. STAR 
measurements [1C§ for <fi px spectra in 0-10%, 10-20%, 
10-20%, 20-30%, 30-40%, 40-50% and 50-60% centrality 
Cu+Cu collisions are shown in FigHb. A Levy function 
also fits the pr spectra in Cu+Cu collisions [lCj . For sim- 
ilar Np ar t, 4> spectra in Au+Au and Cu+Cu collisions are 
similar. The parameters of the Levy function in Au+Au 
collision and Cu+Cu collision are also similar for nearly 
identical for participant numbers. Black lines in FigJHjD 
are the predictions from viscous hydrodynamics. Here 
again, except for the very peripheral (50-60%) collisions, 
data, up to pr=3 GeV are well explained. The results 
indicate that <fi pr spectra, up to pr— 3GeV, both in 
Au+Au and in Cu+Cu collisions are consistent with hy- 
drodynamic evolution of QGP fluid with viscosity over 
entropy ratio 7//s=0.25. The central energy density of 
the fluid in b=0 Au+Au (Cu+Cu) collisions is «5. 1(3.48) 
GeV/fm 3 . 



C. N(n)/N(cj>) vs. p T 

STAR collaboration measured the transverse momen- 
tum dependence of the ratio N (ty / N (<p) in Au+Au colli- 
sions STAR measurements of the ratio in 0-12%, 20- 
40% and 40-60% centrality Au+Au collisions are shown 
in the three panels, a,b,c of Fig[5J The ratio increases 
with pt till pt ~ 3-4 GeV then drops. In peripheral col- 
lisions, the ratio drops at lower pt than in more central 
collisions. Both fl(sss) and 0(ss) are strange particles, 
devoid of any non-strange quarks. The ratio N(fl)/N((f>) 
can shed light on the production mechanism of strange 
particles, specifically, strange baryon and mesons. The 
ratio can also test a model. Correct reproduction of the 
ratio will indicate that the strangeness sector is correctly 
modeled. Models based on recombination of thermal 
strange quarks [llj can reproduce the d> meson pt spec- 
tra up to pt= 5 GeV. The model [1JJ also reproduces 
the ratio N(tt)/N((j>) up to pt ~ 4 GeV, reproduces the 
decreasing trend at pt > 4 GeV. But at pt > 4 GeV, 
the model largely over predict the ratio, indicating that 
at large pt, in recombination models STs are more pro- 
duced than in experiment. In FigO the black lines are 
hydrodynamic predictions for the ratio in Au+Au colli- 
sions. Hydrodynamics predictions are shown up to pt= 
5 GeV. The ratio increases with px, and continue to in- 
crease even at large pt- Evolution of viscous QGP do not 
reproduce the experimental trend that the ratio decreases 
beyond a certain pt- It is not expected also. As noted 
earlier, 4> spectra at pt > 3 GeV are not reproduced. 
In pt range pt < 3 GeV, viscous hydrodynamics appear 
to under predict the ratio in pt range pt < 3 GeV. For 
example in 0-12% centrality collisions, the ratio is under 
predicted by ~ 40%. In more peripheral collisions, the 
ratio is less under predicted. The result is interesting. 
4> pr-spectra, up to pt=3 GeV, are well reproduced in 
viscous hydrodynamics. Apparently, in viscous hydrody- 
namics, fi's are not produced in sufficient number. The 
ratio N(£i)/N(4>) is not measured in Cu+Cu collisions. 
The blue lines in Fig [9] are the predictions for the ratio 
in Cu+Cu collisions. The ratio is nearly identical to that 
in Au+Au collisions. The ratio N(fl)/N((j)) do not show 
any system size dependence. 



D. (f> elliptic flow in Au+Au and Cu+Cu collisions 

STAR collaboration measured 4> meson elliptic flow in 
0-5%, 10-40% and 40-80% and 0-80% (minimum bias) 
centrality Au+Au collisions [8]. STAR measurements 
for the elliptic flow are shown in Fig[l0l In Fig[l0l the 
black lines are the elliptic flow from evolution of viscous 
fluid with central energy density £o=5-l GeV/ fm 3 . Even 
though as mentioned earlier, viscous hydrodynamics is 
not reliable beyond pt—Z GeV {4> pr-spectra are un- 
der predicted), we have shown predictions for flow up 
to pt—5 GeV. In 0-5% collisions, elliptic flow is very 
small, and viscous fluid evolution reproduces the flow. 
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FIG. 9: (color online) Filled circles are STAR data on the pr 
dependence of the ratio N(Q)/N((j>) in Au+Au collisions in 
0-12%, 20-40% and 40-80% centrality Au+Au collisions. The 
black lines are the viscous hydrodynamic predictions for the 
ratio in Au+Au collisions. The blue lines are the predictions 
for the ratio in Cu+Cu collisions. 
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FIG. 10: (color online) Filled circles in panels a,b,c and d 
are STAR measurements for 4> meson elliptic flow in 0-10%, 
10-40% and 40-80% and 0-80% centrality Au+Au collisions. 
The black lines in the figure are the hydrodynamic model 
predictions for elliptic flow in Au+Au collisions. Blue lines 
are the predicted flow in Cu+Cu collisions. 



Elliptic flow in 10-40% collisions is also reproduced in 
the model. However, flow is largely over predicted in 
40-80% centrality collisions. 40-80% centrality collisions 
approximately corresponds to b=ll fm Au+Au collision. 
Hydrodynamic models are not reliable at such peripheral 
collisions. Interestingly, elliptic flow in 0-80% centrality 
collisions is also well reproduced in viscous hydrodynam- 
ics. Considering that only the central energy density 
is fixed to reproduce (j> multiplicity in 0-5% centrality 
Au+Au collisions, reproduction of (j) px spectra and el- 
liptic flow in central and mid central collisions can be 
considered as a great success of viscous hydrodynamics. 



Blue lines in FigfTUl are the predictions for flow in 0- 
5%, 10-40% and 40-80% and minimum bias Cu+Cu col- 
lisions. Compared to Au+Au collisions, elliptic flow is ~ 
10% less in Cu+Cu collisions. It is consistent with our 
predictions for integrated v%. Integrated is also less 
in Cu+Cu than in Au+Au collisions (see FigJTJ) . </> me- 
son elliptic flow in Cu+Cu collisions are not measured yet 
and the predictions can not be tested against experiment. 
Future experiments can verify the predictions. 

Present analysis indicate that QGP fluid, with viscos- 
ity over entropy ratio ry/s=0.25, is consistent with most 
of the published STAR data on <f> production in central 
and mid-central Au+Au and Cu+Cu collisions. The ini- 
tial central energy density of the fluid is «5.1 GeV/ fm 3 
in Au+Au collisions and ps 3.48 GeV/fm 3 in Cu+Cu 
collisions. Present estimate of viscosity and initial en- 
ergy density are obtained by fitting experimental data in 
a hydrodynamic model. Limitations of the model must 
be discussed. We have neglected bulk viscosity. In gen- 
eral bulk viscosity is much smaller than shear viscosity. 
However, recent lattice simulation [|| indicate that trace 
anomaly is non-zero near the cross-over temperature. Us- 
ing the lattice data, Kharzeev et al [4l[ computed bulk 
viscosity of QGP. Near the cross over temperature, bulk 
viscosity can be significantly large. Experimental data 
include effect of both the shear and bulk viscosity. Ne- 
glecting bulk viscosity will result into overestimating the 
shear viscosity. Then ry/s=0.25 is an upper bound of 
QGP viscosity. However, since bulk viscosity is appre- 
ciable only near the cross over temperature, we do not 
expect substantial entropy production due to bulk vis- 
cosity and estimate of shear viscosity will largely remain 
unaltered. 



IV. SUMMARY AND CONCLUSIONS 

To summarise, in the Israel-Stewart's theory of dissi- 
pative hydrodynamics, we have simulated <f> production 
from Au+Au and Cu+Cu collisions at ^=200 GeV. In 
an earlier publication [36j , we have shown that the STAR 
data on <j) mean px in Au+Au collisions is sensitive to 
viscosity and estimated QGP viscosity as 7y/s=0.25. For 
Tj/s=0.25, QGP fluid, thermalised at -^=0.2 fm and ini- 
tialised with central energy density 5.1 GeV/ fm 3 , ex- 
plain centrality dependence of (j) mean px, multiplicity, 
integrated v% and px spectra (up to px ~ 3 GeV). It is 
now shown that the STAR data on cf> elliptic flow in cen- 
tral and mid-central Au+Au collisions are also explained 
in evolution of QGP fluid with viscosity to entropy ratio 
7//s=0.25. rj/s=0.25 is also consistent with STAR data 
on (p meson px spectra in Cu+Cu collisions. In Cu+Cu 
collisions, central energy density is 3.48 GeV / fm 3 , ~1.5 
times less that in Au+Au collisions. We have given pre- 
dictions for elliptic flow in Cu+Cu collisions. Predicted 
elliptic flow in Cu+Cu collisions is ~10% less than that 
in Au+Au collisions. In conclusion, STAR data on in- 
termediate range px < 3 GeV, 4> meson production is 
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Au+Au and Cu+Cu collisions is consistent with hydro- ratio T]/s «0.25. 
dynamic evolution of QGP fluid with viscosity to entropy 
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